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Abstract 



We analyze the impact of the data on long baseline v„ disappearance from the K2K experi- 
ment and reactor v e disappearance from the KamLAND experiment on the determination of the 
leptonic three-generation mixing parameters. Performing an up-to-date global analysis of solar, 
atmospheric, reactor and long baseline neutrino data in the context of three-neutrino oscillations, 
we determine the presently allowed ranges of masses and mixing and we consistently derive the 
allowed magnitude of the elements of the leptonic mixing matrix. We also quantify the maxi- 
mum allowed contribution of Am^ oscillations to CP-odd and CP-even observables at future long 
baseline experiments. 
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I. INTRODUCTION 



Neutrino oscillations are entering in a new era in which the observations from underground 
experiments obtained with neutrino beams provided to us by Nature - either from the Sun 
or from the interactions of cosmic rays in the upper atmosphere- are being confirmed by 
experiments using "man-made" neutrinos from accelerators and nuclear reactors. 



Super-Kamiokande (SK) high statistics data jjj, |2j clearly established that the observed 
deficit in the /i-like atmospheric events is due to the neutrinos arriving in the detector at 
large zenith angles, strongly suggestive of the oscillation hypothesis. This evidence was 
also confirmed by other atmospheric experiments such as MACRO [3] and Soudan 2 0. 
Similarly, the SNO results 0] in combination with the SK data on the zenith angle depen- 

0, SAGE 0, 



dence and recoil energy spectrum of solar neutrinos |6j] and the Homestake 



dence and recoil energy 
GALLEX+GNO 0,E3 



and Kamiokande 



111 ] experiments, put on a firm observational basis 
the long-standing problem of solar neutrinos [12j, establishing the need for v e conversions. 

The KEK to Kamioka long-baseline neutrino oscillation experiment (K2K) uses an 
accelerator-produced neutrino beam mostly consisting of with a mean energy of 1.3 GeV 
and a neutrino flight distance of 250 km toprobe the same oscillations that were explored 
with atmospheric neutrinos. Their results [13] show that both the number of observed neu- 
trino events and the observed energy spectrum are consistent with neutrino oscillations with 
oscillation parameters consistent with the ones suggested by atmospheric neutrinos. 

The KamLAND experiment measures the flux of £/ e 's from nuclear reactors with a energy 
of ~ MeV located at typical distance of ~ 180 km with the aim of exploring with a terrestrial 
beam the region of neutrino parameters that is relevant for the oscillation interpretation 
of the solar data. Their first published 3| results show that both the total number of 
events and their energy spectrum can be better interpreted in terms of v f oscillations with 



14, 



15, 



16, 



12|. 



parameters consistent with the the LMA solar neutrino solution 

Altogether, the data from solar and atmospheric neutrino experiments, and the first 
results from KamLAND and K2K constitute the only solid present-day evidence for physics 
beyond the Standard Model 18j. The minimum joint description of these data requires 
neutrino mixing among all three known neutrinos and it determines the structure of the 
lepton mixing matrix which can be parametrized as [20] 
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(1) 



where = cos 6 \j and sy = sin 6^-. In addition to the Dirac-type phase, 5, analogous to 



that of the quark sector, there are two physical phases associated to the Majorana character 
of neutrinos, which however are not relevant for neutrino oscillations [sjj and will be set to 
zero in what follows. 

In this paper we present the result of the global analysis of solar, atmospheric, reactor 
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and long-baseline neutrino data in the context of three-neutrino oscillations with the aim of 
determining in a consistent way our present knowledge of the leptonic mixing matrix and 
the neutrino mass differences. We make particular emphasis on the impact of the first data 
from long baseline disappearance from the K2K experiment and reactor v e disappearance 
from the KamLAND experiment. 

The outline of the paper is as follows. In Sec. |H] we describe the data included in the 
analysis and we briefly describe the relevant formalism. Sec. 1111 Al contains the results of the 
analysis of the K2K data and their effect on the determination of the parameters associated 
with atmospheric oscillations. We find that the main impact of K2K when combined with at- 
mospheric neutrino data is to reduce the allowed range of the corresponding mass difference. 
When combined with the data from the CHOOZ [22] experiment in a three-neutrino analy- 
sis, this results into a slight tightening of the derived bound on #13 at high CL. In Sec. IIII Bl 
we describe the results from the global analysis including also solar and KamLAND data 
and in Sec. 1111 CI we describe our procedure to consistently derive the allowed magnitude of 
the elements of the leptonic mixing matrix. As outcome of this analysis we also quantify the 
maximum allowed contribution of Am^ oscillations to CP-odd and CP-even observables at 
future long baseline experiments in Sec. IIII Dl . Conclusions are given in Sec. IIVI We also 
present an appendix with the details of our analysis of the K2K data. 



II. DATA INPUTS AND FORMALISM 



We include in our statistical analysis the data from solar, atmospheric and K2K acceler- 
ator neutrinos and from the CHOOZ and KamLAND reactor antineutrinos. 

In the analysis of K2K we include the data on the normalization and shape of the spec- 
trum of single-ring /i-like events as a function of the reconstructed neutrino energy. The total 
sample corresponds to 29 events. In the absence of oscillations 44 events were expected. We 
bin the data in five 0.5 GeV bins with < E rec < 2.5 plus one bin containing all events 
above 2.5 GeV. For QE events the reconstructed neutrino energy is well distributed around 
the true neutrino energy, However experimental energy and angular resolution and more 
importantly the contamination from non-QE events, result into important deviations of the 
reconstructed neutrino energy from the true neutrino energy which we carefully account for. 
We include the systematic uncertainties associated with the determination of the neutrino 
energy spectrum in the near detector, the model dependence of the amount of nQE con- 
tamination, the near/far extrapolation and the overall flux normalization. Details of this 
analysis are presented in the appendix. 

For atmospheric neutrinos we include in our analysis all the contained events from the 
latest 1489 SK data set 1], as well as the upward-going neutrino-induced muon fluxes from 
both SK and the MACRO detector j^]. This amounts for a total of 65 data points. More 
technical description of our simulations and statistical analysis can be found in Refs. 0, E^J ■ 
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We refine our previous analysis j3 of the CHOOZ reactor data and include 
here their energy binned data instead of their total rate only. This corresponds to 14 data 
points (7-bin positron spectra from both reactors, Table 4 in Ref. |22j) with one constrained 
normalization parameter and including all the systematic uncertainties there described. 

For the solar neutrino analysis, we use 80 data points. We include the two measured 
radiochemical rates, from the chlorine and the gallium g Q experiments, the 44 
zenith-spectral energy bins of the electron neutrino scattering signal measured by the SK 
collaboration [6], and the 34 day-night spectral energy bins measured with the SNO 0] 
detector. We take account of the BP00 [26j] predicted fluxes and uncertainties for all solar 
neutrino sources except for 8 B neutrinos. We treat the total 8 B solar neutrino flux as 
a free parameter to be determined by experiment and to be compared with solar model 
predictions. For KamLAND we include information on the observed antineutrino spectrum 
which accounts for a total of 13 data points. Details of our calculations and statistical 
treatment of solar and KamLAND data can be fonnd in Rets. 00. 

In general the parameter set relevant for the joint study of these neutrino data in the 
framework of three-v mixing is six-dimensional: two mass differences, three mixing angles 
and one CP phase. 

Results from the analysis of solar plus KamLAND, and atmospheric data in the framework 
of oscillations between two neutrino states 000,00 imply that the required mass 
differences satisfy that 

Am 2 Q « Am ltm- (2) 

In this approximation the angles 6ij in Eq. can be taken without loss of generality to lie 
in the first quadrant, % G [0, vr/2]. There are two possible mass orderings which we chose 
as 

Am 2 21 = Am 2 G < A «4 - Arr& = > 0; (3) 

Am\ x = Ami < ~ Am li - ~ Am ^ = \ A ™ltJ > 0. (4) 

As it is customary we refer to the first option, Eq. Q, as the normal scheme, and to the 
second one, Eq. (jlj), as the inverted scheme. 

For solar neutrinos and for antineutrinos in KamLAND the oscillations with Am\ 2 ~ 
Amjj are averaged out. The relevant survival probability takes the form: 

Pi: = sin 4 9 13 + cos 4 ^C(AmL M , (5) 

where P e 2 e 1/ (Am| 1 , f9 12 ) is the survival probability for 2v mixing, which, for solar neutrinos, 
is obtained with the modified sun density N e — > cos 2 9i^N e . So the analysis of solar and 
KamLAND data depends on three of the five oscillation parameters: Am| 1 ,f9i2 and #13. 

Conversely for small Am^ the three-neutrino oscillation analysis of the atmospheric 
and K2K neutrino data can be performed in the one mass scale dominance approximation 
neglecting the effect of Am^. In this approximation the angle 6 , 12 can be rotated away 
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and it follows that the atmospheric and K2K data analysis restricts three of the oscillation 
parameters, namely, Am^ = Am^, # 23 and # 13 and the CP phase 5 becomes unobservable. 
The Vu survival probability at K2K is 



pK2K _ I _ a ( 4 2 2 , 2 2 2 \ - 2 ( Am 32 r \ 
ju/i _ 1 4 l S 23 S 13 C 13 + C 13 S 23 C 23j Sm I |^ ^ J 



2 COS 26*23 . /-, 2 cos 2^23 n n K2K,2vf a 2 /. \ /-/V 4 \ / c \ 

- s i3— 5 + ( 1 - S i3— 72 ) P w (Am 32 ,^ 3 2) + 0(s 13 ) . (6) 

C 23 C 23 

For atmospheric neutrinos in the general case of three-neutrino scenario with # 13 7^ 
the presence of the matter potentials become relevant. We solve numerically the evolution 
equations in order to obtain the oscillation probabilities for both e and \x flavours, which are 
different for neutrinos and anti-neutrinos. Because of the matter effects, they also depend 
on the mass ordering being normal or inverted. In our calculations, we use for the matter 
density profile of the Earth the approximate analytic parametrization given in Ref. of 
the PREM model of the Earth [29]. 

he reactor neutrino data from CHOOZ provides information on the survival probabil- 

i3. 



ity 



24, 



30 
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The second equality holds in the approximation Am^ <C E„/L which can be safely made 
for the presently allowed values of Am| x (rSI Thus the analysis of the CHOOZ reactor 
data involves only two parameters: Am 32 and the mixing angle #13. 

In summary, oscillations at solar+KamLAND on one side, and atmospheric+K2K oscil- 
lations on the other side, decouple in the limit #13 = 0. In this case the values of allowed 
parameters can be obtained directly from the results of the analysis in terms of two-neutrino 
oscillations and normal and inverted hierarchies are equivalent. Deviations from the two- 
neutrino scenario are determined by the size of the mixing # 13 . 

The allowed ranges of masses and mixing obtained in our two-neutrino oscillation analysis 
of solar+KamLAND data can be found in Ref. ^| and we do not reproduce them here. We 
discuss next the results of our analysis of K2K data and its impact on the determination of 
the parameters relevant in atmospheric oscillations. 



III. RESULTS 

A. Am 32 Oscillations: Impact of K2K Data 

For sake of comparison with the K2K oscillation analysis, we discuss first the results of 
our analysis of K2K data for pure — ► v T oscillations which are graphically displayed in 
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Fig. Ill 

We show in the left panel of Fig. ^ the allowed region of (Am 2 , sin 2 (28)) from our analysis 
of K2K data. The best fit point for this analysis is at Am 2 = 2.7 x l(T 3 eV 2 , sin 2 (26 l ) = 0.92 
with xLin = 9-3 (the corresponding best fit as obtained from the K2K collaboration is a 
Am 2 = 2.7 x 10 _3 eV 2 , sin 2 (2#) = 1). We notice that the non-maximality of the mixing 
angle in our analysis is not statistically significant as maximal mixing is only at A% 2 = 
0.15. The energy spectrum for this point is shown in the right panel together with the 
data points and the expectations in the absence of oscillations. Our results show very 
good agreement with those obtained by the K2K collaboration [l3j. Also displayed in the 
figure are the corresponding regions from our latest atmospheric neutrino analysis j^j]. 
As seen in the figure, the K2K results confirm the presence of neutrino oscillations with 
oscillation parameters consistent with the ones obtained from atmospheric neutrinos studies. 
Furthermore, already at this first stage, it provides a restriction on the allowed range of Am 2 , 
while their dependence on the mixing angle is considerably weaker. 




FIG. 1: — > v T oscillation analysis of K2K data. On the left panel we show the allowed (2d.o.f) 
regions on Am 2 , sin 2 29 at 90% CL (solid) and 99%CL (dotted). The best fit point is marked with 
a thick dot. The shadowed regions are the 90, and 99% CL of the atmospheric neutrino analysis 
with best point marked by a star. The right panel shows the spectrum of K2K events as a function 
of the reconstructed neutrino energy for the 6 bins used in the analysis. The data points are shown 
together with their statistical errors. The dotted histogram and the shadow boxes represent our 
prediction and the K2K MC prediction in the absence of oscillations respectively. The full line 
represents the expected distribution for the best fit point for — > v T oscillations. 

In the framework of 3u mixing the analysis of K2K, atmospheric and CHOOZ data 
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FIG. 2: 3u oscillation analysis of the atmospheric, CHOOZ and K2K data. The left, center and right 
panels show the dependence of A% 2 on Am| 2 , tan 2 #23 an d sin 2 #13 for the analysis of atmospheric, 
CHOOZ and K2K data (full line) compared to the previous bound before the inclusion of K2K 
(dotted line). The individual 1 (3) o bounds in Eq. @ can be read from the figure with the 
condition A% 2 < 1 (9). 

provides information on the parameters Am^, #23, and #13. We define: 

XATM+CHOOZ+K2K(A m 32> 6*23, #13) = XATM( Am 32, #23, #13) 

+XcHOOz( Am 32> M + x! 2 K( Am L #23, M • (8) 

In the three panels of Fig. El we show the bounds on each of the three parameters obtained 
from this analysis (full lines). For comparison we also show the corresponding ranges for 
the analysis of atmospheric and CHOOZ data alone (dotted lines). The corresponding 
subtracted minima are given in Table |1J The results in the figure are shown for the normal 
mass ordering, but once the constraint on # 13 from CHOOZ is included in the analysis, the 
differences between the results for normal and inverted mass ordering are minimal. The 
careful reader may notice that the \ 2 P er d.o.f. seems too good. As seen in Table [Q this 
effect is driven by the atmospheric data and it was already the case for the previous SK data 
sample. It is partly due to the very good agreement of the multi-GeV electron distributions 
with their no-oscillation expectations. However, as discussed in Ref. J^, Xmm 1S om Y 
below its characteristic value, not low enough to be statistically suspected. 

In each panel the displayed \ 2 has been marginalized with respect to the other two 
parameters. From the figure we see that the inclusion of K2K data in the analysis results 
into a reduction of the allowed range of Am 2 2 while the allowed range of #23 is not modified. 
The reduction is more significant for the upper bound of Am 2 2 while the lower bound is 
slightly increased. More quantitatively we find that the following ranges of parameters are 
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ATMOS ATMOS+CHOOZ ATM0S+CH00Z+K2K 


Data points 

Am§ 2 
eV 2 

tan 2 2 3 
sin 2 0i3 


65 65+14=79 65+14+6=85 
2.7 x 1(T 3 2.55 x 1(T 3 2.6 x 1(T 3 
1 1 1 
0.015 0.009 0.009 
39.7 45.8 55.1 



TABLE I: Minimum \ 2 values and best-fit points for the 3u oscillation analysis of atmospheric, 
CHOOZ and K2K data. 

allowed at la (3a) CL from this analysis 

(1.5)2.2 < Amyi(T 3 eV 2 < 3.0(3.9), 

(0.45)0.75 < tan 2 ^3 < 1.3(2.3). (9) 

These ranges are consistent with the results from the 2-neutrino oscillation analysis of K2K 
and atmospheric data in Ref. j^ . 

Concerning the "generic" 3z/ mixing parameter # 13 , Eq. © shows that its effect on both 
the normalization and the shape of the spectrum is further suppressed near maximal 
mixing by cos 2023 ~ 0. As a consequence K2K alone does not provide any bound on 
#13. However Fig. |21c illustrates how the inclusion of the long-baseline data results into 
a tightening of the bound on # 13 (at large CL) when combined with the atmospheric and 
CHOOZ data. This is an indirect effect due to the increase in the lower bound on Am| 2 ■ 
In the favoured range of Am 2 2 the oscillating phase at CHOOZ is small enough so that it 
can be expanded and the oscillation probability of v e depends quadratically on Am^. As a 
consequence the bound on the mixing angle from CHOOZ is a very sensitive function of the 
allowed values for Am| 2 . The increase of the lower bound on Am 2 2 due to the inclusion of 
the K2K data leads to the tightening of the derived limit on #13 at high CL. From Fig. [21c 
and Table |U we also see that the best fit point is not exactly at sin 2 #13 = 0, although this 
is not very statistically significant. This effect is due to the atmospheric neutrino data. In 
particular to the slight excess of sub GeV e-like events which is better described with a 
non- vanishing value of #13. 



B. Global Analysis 

We calculate the global \ 2 by fitting all the available data 

Xglobal( A ™21, Am 32, 012, 023, #13 ) = Xsolar^^l, #12, M + XKland( A ™21, #12, M 

+Xatm( Am 32, #23, #13 ) + xLk^^L #23, #13 ) 
+XCHOOZ( A ^32 ? 013) • (10) 



8 



The results of the global combined analysis are summarized in Fig. |3] and Fig. |U in which 
we show different projections of the allowed 5-dimensional parameter space. 

In Fig. El we plot the individual bounds on each of the five parameters derived from the 
global analysis (full line). To illustrate the impact of the K2K and KamLAND data we also 
show the corresponding bounds when K2K is not included in the analysis (dotted line) and 
when KamLAND is not included (dashed line). In each panel the displayed x 2 has been 
marginalized with respect to the other four parameters. The subtracted minima for each of 
the curves are given in Table ITT1 

Figure 01 illustrates that the dominant effects of including KamLAND are those derived in 
the two-neutrino oscillation analysis of solar and KamLAND data 

0,1 Qjfl: the determination 

of Am|i (panel a) to be in the LMA region and a very mild improvement of the allowed 
mixing angle #12 (panel b). In other words, the inclusion of the 3z/ mixing structure in the 
analysis of solar and KamLAND data does not affect the determination of these parameters 
once the additional angle 6*i 3 is bounded to be small. The slight tightening of the #13 limit 
due to the inclusion of K2K data does not have any impact in the determination of the 
bounds on Am^j and $12- Quantitatively we find that the following ranges of parameters 
are allowed at la (3a) CL from this analysis 

(5.4)6.7 < Am^/lO^eV 2 < 7.7(10.) and (14.) < Am^/lO^eV 2 < (19.), 

(0.29)0.39 < tan 2 ^i 2 < 0.51(0.82) (11) 

The range of Am^ on the right of the first line in Eq. (jllj) correspond to solutions in the 
upper LMA island (see Fig. EJa). At present the results of the solar and Kamland analysis 
still allow for this ambiguity in the determination of Am^ at CL > 2.5a. This reflects the 
departure from the parabolic (gaussian) behaviour of the Am^ dependence of Xgiobai and 
the presence of a second local minima. With improved statistics KamLAND will be able to 
resolve this ambiguity 17, 3^|. 



Comparing the full line on the #13 panel in Fig.Elc with the corresponding one in Fig. EJc 
we see that the inclusion of the solar+KamLAND data does have an impact on the allowed 
range of #13. However the comparison of the full and dashed line in Fig. Elc illustrates that 
the impact is due to the solar data. Eq. (jSJ) shows that a small # 13 does not significantly 
affect the shape of the measured spectrum at KamLAND. On the other hand, the overall 
normalization is scaled by cos 4 #13 and this factor has the potential to introduce a non- 
negligible effect (in particular in the determination of the mixing angle 612 (3)- Within its 
present accuracy, however, the KamLAND experiment cannot provide any further significant 
constraint on #13. Altogether, the derived bounds on #13 from the global analysis are: 

sin 2 #13 < 0.02 (0.052) , (12) 

at la (3a). 

Finally, comparing Fig. Old and Fig. |31e with the corresponding curves in Fig. [21a and 
Fig. [2]b, we see that the additional restriction on the possible range of #13 imposed by 
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FIG. 3: Global 3v oscillation analysis. Each panels shows the dependence of Ax 2 on each of the 
five parameters from the global analysis (full line) compared to the bound prior to the inclusion of 
the K2K (dotted line) and KamLAND data (dashed line). The individual 1 (3) a bounds in Eqs. 
(|TT|) and (|T2|) @ can be read from the corresponding panel with the condition A% 2 < 1 (9). 

the solar data does not quantitatively affect the dominant effect of the inclusion K2K -the 
improved determination of Am| 2 - Thus the allowed ranges of Am§ 2 and # 23 in Eq. © are 
valid for the global analysis as well. 

The ranges in Eqs. 0, (fTTj) and (fT2"j) are not independent. In Fig.0]we plot the correlated 
bounds from the global analysis for each pair of parameters. The regions in each panel are 
obtained after marginalization of Xgiobai i n Eq (flOf with respect to the three undisplayed 
parameters. The different contours correspond to regions defined at 90, 95, 99 % and 3<r 
CL for 2d.o.f (A% 2 = 4.61,5.99,9.21,11.83) respectively. From the fi gure we see that the 
stronger correlation appears between 9 13 and Ar??| 2 as a reflection of the CHOOZ bound. 
In general, because of the correlations, the ranges in Eqs. 0, (fTTj) and (fT2"j) cannot be 
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tarftf 



23 



0.05 0.1 

sin 2/ $ 13 



FIG. 4: Global 2>u oscillation analysis. Each panels shows 2-dimensional projection of the allowed 
5-dimensional region after marginalization with respect to the three undisplayed parameters. The 
different contours correspond to the two-dimensional allowed regions at 90, 95, 99 % and 3cr CL. 



directly used in deriving the corresponding entries in the U mixing matrix as we discuss 
next. 
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GLOBAL 


GLOBAL-K2K GLOBAL-KLAND 


Data points 


178 


172 


165 


Am 21 
eV 2 

tan 2 #12 

Am§ 2 
eV 2 


7.1 x 1(T 5 

0.45 
2.6 x 1(T 3 


7.1 x 10- 5 

0.45 
2.5 x 10~ 3 


5.8 x 10- 5 

0.45 
2.6 x 10" 3 


tan 2 6*23 


1 


1 


1 


sin 2 6<i3 


0.009 


0.009 


0.009 


x 2 - 

A,min 


136 


127 


130 



TABLE II: Minimum x 2 values and best-fit points for the global 3v oscillation analysis. 



C. Determination of the Leptonic Mixing Matrix 



We describe in this section our procedure to consistently derive the allowed ranges for 
the magnitude of the entries of the leptonic mixing matrix. We start by defining the mass- 
marginalized x 2 function: 

Xmix,global(#12, #23, #13) = min^^ ^m^Xglobal (^ m 21 > ^ m 32> #12, #23, #13) • (13) 

We study the variation of Xmix.giobai as function of each of the mixing combinations in U 
as follows. For a given magnitude C/^ of the entry U(i,j) we define x 2 {Uij) as the mini- 
mum value of Xmix, g iobai(#i2, #23, #13) with the condition \U(i, j)(#i2, #23, #13)! = U itj . In this 
procedure the phase 5 is allowed to vary freely between and n. The allowed range of the 
magnitude of the entry ij at a given CL is then defined as the values U it j verifying 

XW) - Xglobal.min < A* 2 (CL, Id.O.f) , (14) 

with Xgiobai.min = 136- This is equivalent to having done the full analysis in terms of the 
independent matrix elements - of which, in the hierarchical approximation, only three are 
experimentally accessible at present (and can be chosen for instance to be \U e 2\, \U e s\, and 
|£^3|)~ and find the allowed magnitude of each \Uij\ by marginalization of 

Xmix,global(l^e2|, \U e ?\, \U ^\) , (15) 

with the use of unitarity relations and allowing a free relative phase 5. 

With this procedure we derive the following 90% (3a) CL limits on the magnitude of the 
elements of the complete matrix 



U = 



/(0. 73) 0.79 to 0.86 (0.88) (0.47) 0.50 to 0.61 (0.67) to 0.16 (0.23) \ 
(0.17) 0.24 to 0.52 (0.57) (0.37) 0.44 to 0.69 (0.73) (0.56) 0.63 to 0.79 (0.84) 
v (0.20) 0.26 to 0.52 (0.58) (0.40) 0.47 to 0.71 (0.75) (0.54) 0.60 to 0.77 (0.82) j 



(16) 
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By construction the derived limits in Eq. (|16p are obtained under the assumption of the 
matrix U being unitary. In other words, the ranges in the different entries of the matrix 
are correlated due to the fact that, in general, the result of a given experiment restricts 
a combination of several entries of the matrix, as well as to the constraints imposed by 
unitarity. As a consequence choosing a specific value for one element further restricts the 
range of the others. 



D. Am^ oscillations at future LBL experiments 

In general, correlations between the allowed ranges of the parameters have to be consid- 
ered when deriving the present bounds for any quantity involving two of more parameters. 
This is the case, for example, when predicting the allowed range of CP violation at future 
experiments as discussed in Ref. |3o1 |. 

Here we explore the possible size of effects associated with Am^ oscillations (both CP 
violating and CP conserving) at future LBL experiments to be performed either with con- 
ventional superbeams (conventional meaning from the decay of pions generated from a 
proton beam dump) at a Neutrino Factory with neutrino beams from muon decay in 
muon storage rings . 

The "golden" channel at these facilities involve the observation of either "wrong-sign" 
muons due to v e — > (or u e — > z/ M ) oscillations at a neutrino factory or or the detection of 
electrons (positrons) due to — > v e {p^ — > u e ) at conventional superbeams . In either case 



the relevant oscillation probabilities in vacuum are accurately given by 38|, 3^ 



F„ e „„ = 4 sin 2 2^ sin 2 [~^-J + 4 sin 2 2# 12 [~^- 

+ Jcos (, + Si „ (^M) = P- + i~ + , (17) 

with J = C13 sin 2# 12 sin 26 13 sin 2# 2 3- P so1 contains the contribution to the probability due 
to longer wavelength oscillations while P inter gives the interference between the longer and 
shorter wavelength oscillations and contains the information on the CP-violating phase 5. 
In order to quantify the present bounds on these contributions we factorize the baseline and 
energy independent parts as: 

psol = (F sol )2 ( _M 2 pinter = pinter f & _ { ( ( ± 

v ' \AE J V 4:E J \ 4E J \AE / 

E so1 = c 23 sin 2#i 2 Am 2 21 E intei = c 13 sin 2# 12 sin 29 13 sin 2# 23 Am 2 21 . (18) 

For very long baselines, for which the presence of matter cannot be neglected, the expressions 
above for F so1 and ,F mter still hold as the coefficients of the dominant contributions to the 
probabilities in the expansion in the small parameters #13 and Am^ 
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5 1015 20 1 2 3 4 5 



F sol ( 1 -5 eV 2 ) Am 2 32 ( 1 -3 eV 2 ) 

FIG. 5: (a) Dependence of A% 2 on F so1 . (b) Allowed regions of F inter versus Am| 2 at 90, 95, 99% 
and 3<T. The best value is marked with a star. 

We show in Fig. EJthe present bounds on the coefficients F so1 and _p mter . In general the 
dependence on Am^ of the interference term cannot be factorized because, depending on the 
considered baseline and energy, the oscillating phase with Am§ 2 may be not small enough 
to be expanded. For this reason we show in Fig. 0b the 2-dimensional allowed region of 
pinter versus Am 2 2 . In the figure we mark with a star the best value for F mter as obtained 
from this analysis which is not vanishing due to the small but non-zero best fit value of 
sin #13. This is, however, not statistically significant as F mter = is at A% 2 = 0.9. The 
negative slope in the upper part of the 90% and 95% CL region in Fig. EDb is a reflection of 
the anti-correlation between the Am^ and sin 2 # 13 constraints from the CHOOZ experiment 
(see Fig. Bf) 

From this study we find the following la (3a) [ld.o.f] bounds 

F sol /(10- 5 eV 2 ) = 4.6 ± 0.6 (4.6± 2 ;^ and 9.5l|^) , (19) 
< F intcr /(10- 5 eV 2 ) < 1.9 (5.5) , (20) 

where the bounds on .F mter are shown for the best fit value of Am 2 2 = 0.0026 eV 2 . The larger 
values for the 3a range in Eq. (fT9"|) and _F mter correspond to solutions of the solar+KamLAND 
analysis lying in the higher Am^ island (see Fig. EJa and discussion below Eq. (fTTf ). 

IV. SUMMARY AND CONCLUSIONS 

We have presented the results of an updated global analysis of solar, atmospheric, reactor 
and long-baseline neutrino data in the context of three-neutrino oscillations making special 
emphasis on the impact of the recent long baseline v ^ disappearance data from the K2K 
experiment and reactor v e disappearance from the KamLAND experiment. We find that 
the dominant effect of the inclusion of K2K and KamLAND data is the reduction of the 
allowed range of Am^ and Am 21 respectively while the impact on the mixing angles #23 
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and 612 is marginal. The increase of the lower bound on Am| 2 due to the inclusion of the 
K2K data leads also to a slight tightening of the derived limit on # 13 at high CL. Our results 
on the individual allowed ranges for the oscillation parameters are given in Eqs. (jHJ), 
and (fT2"|) and graphically displayed in Fig. El The correlations between the derived bounds 
are illustrated i in Fig. 0J As outcome of the analysis we have presented in Eq. (jl6j) our 
up-to-date best determination of the magnitude of the elements of the complete leptonic 
mixing matrix. Finally we have quantified the allowed contribution of Am^ oscillations to 
CP-odd and CP-even observables at future long baseline experiments with results presented 
in Fig. El and Eqs. and (J271j) . 



APPENDIX A: ANALYSIS OF K2K DATA 

In this appendix we describe our calculation of the K2K spectrum and our statistical 
analysis of the K2K data jlil ]. 

We use in our statistical analysis the K2K data on the spectrum of single-ring /i-like 
events. K2K present their results as the number of observed events as a function of the 
reconstructed neutrino energy. The reconstructed neutrino energy is determined from the 
observed /1 energy in the event, E^, and its scattering angle with respect to the incoming 
beam direction, cos 9^, as 

E,« = , (21) 

m N - E^+ cos dp 

where rriN is the nucleon mass. In Fig. Q]we show their data binned in five 0.5GeV bins 
with < E TCC < 2.5 plus one bin containing all events above 2.5 GeV. The total sample 
corresponds to 29 events. In the absence of oscillations 44 events were expected. 

For QE events, v^n — > fip , and assuming perfect and cos 9^ determination, E rec = E v . 
Experimental energy and angular resolution, nuclear effects, and, more importantly, the 
contamination from non-QE (nQE) events, v^N — > fiX, in the sample, result into important 
deviations of the defined E rec from the real E v . From simple kinematics one finds that in 
nQE events there is a shift in the reconstructed neutrino energy with respect to the true 
neutrino energy 



E = E 

-^rec 



M x -m% 



-1 



<E„, (22) 

m N E^ - _ 

where Mx is the invariant mass of the hadronic system produced together with the muon in 
the interaction. At the K2K energies the most important nQE contamination comes from 
single pion production, which occurs via the A resonance. At the largest energies there is a 
small contribution from deep inelastic scattering. 

Thus in general the observed spectrum of single-ring /i-like events in K2K can be obtained 
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as 



N th (E rcc ) = N norm J ^ SK {E v )P^{E v )[a QE {E v )e^{E 1/ )r QE {E v ,E Tec ) + 



fnQE <y n QE{E v ) e n Q E (E u ) r n Q E (E u , E TCC )~\dE u , (23) 

where $sk(E v ) is the expected spectrum at the SK site in the absence of oscillations. 
PpupiEv) is the survival probability of z/ M for a given set of oscillation parameters. f n QE = 0.93 
is a rescale factor of the expected contamination from nQE events as obtained from MC sim- 
ulation by the K2K collaboration 13]. <jQE{nQE){E u ) are the neutrino interaction cross sec- 

1-/? 

tions. Zqeuqe) are ^ ne detection efficiencies for 1-ring /x-like events at SK. tqe{uQE){E u , E Tec ) 
are the functions relating the reconstructed energy and the true neutrino energy. N norm is the 
normalization factor which is chosen so that in the absence of oscillations the total integral 
gives 44 events. 

In our calculation we use the neutrino spectrum &sk as provided by the K2K collab- 
oration j^J, 41 1 . This flux was estimated from the flux measured in the near detector by 



multiplying it by a MC simulated ratio of the fluxes between the near and far detectors. We 
further assume that the detection efficiencies for 1-ring //-like events at SK are the same for 
the K2K analysis as for the atmospheric neutrino analysis (further details and references 
can be found in Ref. [3]). 

At present there is not enough information from the K2K collaboration on the 
r QE(nQE)(E u , E Tec ) functions. In our calculation we have used a physically-motivated form 
for those functions. We include in the functions ^QE(nQE)(E u , E Tec ) the dominant effect in 
the missreconstruction of the neutrino energy: the shift in the reconstructed neutrino energy 
due to the different kinematics of the nQE events as described by Eq. (J22J)- We also include 
the (subdominant) effects due to the experimental energy and angular resolutions which 
smear the measured muon energy and angle 6^ around their true values E'^ and 9'^. 

r QE {E v , E rcc ) = — -L-j / dE'^dE,de daQE{ ^ E ^ Res E (E, - E'^Res^ - fy8{E v - E' vec ) 
r nQE (E u ,E iec ) = ^— / dE'^dBdM x ^^fS ReSE (^ " K) Res ^ ~ 



S I E u -E> ec 



^ M 2 x -m 2 N n 
m N E>^ - 



(24) 
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where 



EL 



E n 



m N E^ - Si m N -E' +p' cos ^ 



TOjv - Ef, + p fl cos9'^ l 
Resv{K-E v ) = 

Res e (6' u - 9 V ) 



2naE 
1 



(25) 



Following SK data [2] we use an energy resolution for the muons of a E /E^ = 3% and 
an angular resolution ag = 3° (see also Ref. for further details). Notice that in the 
expressions above the true angle of the muon, 6' , is not an independent variable but it is 
related by the kinematics of the process to the initial neutrino energy, E v , the final muon 
energy, E' and the invariant mass of the hadronic system, Mx- The final result on the 



number of expected events in each E res bin is obtained by substituting Eqs. (|2^j) and (J23j) in 
Eq. (J23)) and numerically integrating for the kinematical variables in the corresponding range 
of E ves . In this procedure the only free parameter to adjust is the overall normalization. The 
shape of the spectrum is then fully determined. 

In order to verify the quality of our simulation we compare our predictions for the energy 
distribution of the events with the Monte Carlo simulations of the K2K collaboration in 
absence of oscillation. In Fig. we show our predictions superimposed with those from the 
experimental Monte Carlo (obtained from Fig. 2 in Ref. |l3|), both normalized to the 44 
expected events in the absence of oscillations. The boxes for the MC prediction represent 
the systematic error bands. We can see that the agreement in the shape of the spectrum is 
very good. 

In our statistical analysis of the K2K data we use Poisson statistics as required given the 
small number of events. We include the systematic uncertainties associated with the deter- 
mination of the neutrino energy spectrum in the near detector (ND), the model dependence 
of the amount of nQE contamination parameter f n QE, the near/far extrapolation (F/N) and 



the overall flux normalization (NOR) 4£], |4l| . The errors on the first three items depend 



on energy and have correlations am ong the different energy bins. We account for all these 



effects by using the \ 2 function 32J, |40 , 



41 



XK2K 



miiij 



D 

[ 2 E 



i=l 
7 



N th _ jyexp _ ^yexp 



N, 



j\r ex P 



_|_ jND,nQE^ND,nQE^-l jND,nQE _|_ jNOR 2 



j,k=l 



(26) 



where iVf = iV? h (l + ff /N af /N + f» D af D + f"QE a f E + fNOR a NORy B y^ inf we denote 
the minimization with respect to the systematic shift parameters (or pulls |32j |) f^2^ 6 , f£Lifii 
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f n Q E ; and f NOR We use the systematic errors and their correlations as provided by K2K 



collaboration 



13, 



40 



41 



For instance 



a NOR = 5% , 

af /N = 2.5% , 4.3% , 6.5% , 10.4% 11.1% , 12.2% , 

a ND = 4g% ^ ?1% ; Q% ^ 7 1% g>4% ; n l% ^ (2?) 



a nQE 



13% , 8.9% , 6% , 3.8% 3.% ,5.5 



for z = l,6 respectively. 

Thus in our analysis we use both the shape and the normalization of the 29 single-ring /z- 
like events. In their analysis K2K uses only the spectrum shape (but not the normalization) 
of the 29 single-ring /i-like events plus the overall normalization of their total sample of fully 
contained events (a total of 56). We cannot use the normalization from the additional 27 
events in the lack of more detailed information from the K2K collaboration on the efficiencies 
for multi-ring events. Nevertheless, as described in Sec. IIII Al the results of our oscillation 
analysis are in good agreement with those from the K2K analysis. 
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